Feynman gauge on the lattice: 
new results and perspectives 
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Abstract. We have recently introduced a new implementation of the Feynman gauge on the lattice, 
based on a minimizing functional that extends in a natural way the Landau-gauge case, while 
preserving all the properties of the continuum formulation. The only remaining difficulty with our 
approach is that, using the standard (compact) discretization, the gluon field is bounded, while 
its four-divergence satisfies a Gaussian distribution, i.e. it is unbounded. This can give rise to 
convergence problems when a numerical implementation is attempted. In order to overcome this 
problem, one can use different discretizations for the gluon field, or consider an SU(A^e) group with 
sufficiently large Nc- Here we discuss these two possible solutions. 
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INTRODUCTION 

The behavior of Green's functions in the infrared limit of Yang-Mills theories should 
give us some insights into the low-energy properties of these theories. Since these 
functions depend on the gauge condition, considering different gauges could help us gain 
a better understanding of the (non-perturbative) low-energy hallmarks of QCD, such as 
color confinement. In the last 20 years, several groups have used lattice simulations to 
study propagators and vertices of Yang-Mills theories in Landau gauge [1], Coulomb 
gauge [2, 3, 4], A-gauge (a gauge that interpolates between Landau and Coulomb) [5, 6] 
and maximally Abelian gauge [7, 8]. 

On the other hand, until recently, the numerical gauge fixing for the linear covariant 
gauge — which is a generalization of Landau gauge — was not satisfactory [9, 10, 
11, 12, 13, 14, 15, 16]. In Ref. [17] we have introduced a new implementation of the 
linear covariant gauge on the lattice that solves most problems encountered in earlier 
implementations (see [17, 18] for a short review of early works). As explained in the 
abstract, the only problem still affecting our method, as well as any formulation of 
the linear covariant gauge on the lattice, is due to the fact that the gluon field (x) 
is bounded in the usual compact formulation of lattice Yang-Mills theories. On the 
contrary, the functions A'^(x) satisfy a Gaussian distribution, i.e. they are unbounded. 
Thus, one has to deal with convergence problems [19] when numerically fixing the gauge 



condition 

d^,A^^{x)=A\x). (1) 

Since the real- valued functions A^(x) are generated using a Gaussian distribution with 
width it is clear that this problem becomes more severe when t, is larger and/or 
when the lattice volume is larger. Here we discuss two possible solutions for this 
problem, namely we consider different discretizations for the gluon field or a gauge 
group SU(A^c) with sufficiently large Nc- 

LINEAR COVARIANT GAUGE ON THE LATTICE 

We want to impose the gauge condition (1) on the lattice. Landau gauge, which corre- 
sponds to the case A*(x) = 0, is obtained on the lattice by minimizing the functional 

'^Lcm = -TrY,8{x)U^,{x)g\x + e^,) . (2) 

Here Un{x) are link variables and g{x) are site variables, both belonging to the S\J(Nc) 
group. The sum is taken over all lattice sites x and directions ju. Also, Tr indicates trace 
in color space. For the linear co variant gauge we can look for a minimizing functional 
of the type S'icgW^ ,g,A\. If recalls that solving the system of equations 50 = c, 
where 5 is a matrix and and c are vectors, is equivalent to minimizing the quadratic 
form j^B^ — ^c, then it is obvious that in our case we should look for a minimizing 
functional of the type [17] 

c^LCG[U«,8,A]r^c^LGm-gA. (3) 

Indeed, the lattice linear co variant gauge condition can be obtained by minimizing^ 

<^LCG[U^g, A] = Ag[U'] + 9? rr ig{x) A{x) , (4) 

X 

where S'lg[U^] is defined above in Eq. (2) and 9? indicates real part. This can be checked 
by considering a one-parameter subgroup g{x, t) = exp [/Ty*(x) A^] . Here we indicate 
with A* a basis for the SU(A'c) Lie algebra and with y*(jc) any real- valued functions. 
Indeed, the stationarity condition implies the lattice linear covariant gauge condition^ 

E4w-^m(^-^m) = a^^). (5) 

At 

Also, the second variation (with respect to the parameter t) of the term ig{x)A{x) is 
purely imaginary and it does not contribute to the Faddeev-Popov matrix i.e. ^ 



One should stress that, in the minimization process, the link variables Un (x) get gauge-transformed to 
g{x)Un{x)g\x + en), while the A* (x) functions do not get modified. 

^ Note that periodic boundary conditions yields Ljc^^(^) = 0- This equality has to be enforced explicitly, 
within machine precision, when the functions A*(x) are generated. 



TABLE 1. Smallest value of j3 for 
which the numerical gauge-fixing 
algorithm showed convergence. Re- 
sults are reported for the three dif- 
ferent discretizations and for five 
different values of the gauge param- 
eter ^ . 
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is a discretized version of the usual Faddeev-Popov operator —d -D. Let us note that 
having a minimizing functional for the linear covariant gauge implies that the Faddeev- 
Popov operator ^ is positive-definite and that the set of its local minima defines the 
first Gribov region Q.? 

It is interesting to note that one can interpret the Landau-gauge functional Sic [f/^] as a 
spin-glass Hamiltonian [21] for the spin variables g{x) with a random interaction given 
by U^{x). Then, our new functional corresponds to the same spin-glass Hamiltonian 
when a random external magnetic field A{x) is applied. 

Note also that the functional (f^ccf^'^,,?:^] is linear in the gauge transformation 
{g{x)}. Thus, one can easily extend to the linear covariant gauge the gauge-fixing 
algorithms usually employed in the Landau case [22, 23, 24]. We refer the reader to 
References [17, 18] for tests of convergence of the numerical gauge fixing. There we 
have also checked that the quantity Di{p^)p^, where Di{p^) is the longitudinal gluon 
propagator, is approximately constant for all cases considered, as predicted by Slavnov- 
Taylor identities. This verification failed in previous formulations of the lattice linear 
covariant gauge [13, 16]. 

DISCRETIZATION EFFECTS 

As explained above, the standard discretization of the gluon field (x) is bounded. 

Since the functions A*(x) are generated using a Gaussian distribution, it is clear that 
Eq. (5) cannot be satisfied if A*(x) is too large. A possible solution to this problem is 
to use different discretizations of the gluon field. We did some tests in the SU(2) case 
using the angle (or logarithmic) projection [25] and the stereographic projection [26] 
(for a slightly different implementation of the stereographic projection see also [27]). 
Note that, in the latter case, the gluon field is unbounded even for a finite lattice spacing 
a. 



^ This region has been studied analytically in [20], for a small value of ^ , but a similar numerical study is 
still lacking. 



TABLE 2. Values of the lattice cou- 
pling /3 considered for the gauge 
groups SU(2), SU(3) and SU(4) with a 
gauge parameter <^ = 1 . 



Nc 


/3i 


ft 


ft 


ft 


2 


3.0 


2.485 


2.295 


2.44 


3 


6.75 


6.67 


6.07 


5.99 


4 


12.0 


12.59 


11.43 


10.97 



In particular, considering the standard discretization, the angle projection and 
the stereographic projection for the lattice volume V = 8'*, gauge parameter 
£, = 0.01,0.05,0.1,0.5,1.0 and lattice coupling j8 = 2.2, 2.3,..., 2.9, 3.0 we checked 
(using for the numerical gauge fixing the so-called Cornell method [22, 23, 24]) in 
which cases we were able to fix the covariant gauge condition effectively. Our results, 
reported in Table 1, clearly show that the angle projection is already an improvement 
compared to the standard discretization and that the best convergence is obtained when 
using the stereographic projection. 



CONTINUUM LIMIT 



Note [16] that the continuum relation 



(6) 



can be made dimensionless — working in a generic J-dimensional space — by multi- 
plying both sides by a^go- Since j8 = 2Nc/{a^~^gQ) [in the SU(A^c) case], we have that 
the lattice quantity 



x h '- 
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(7) 



becomes 



1 1 



b 



(8) 



in the formal continuum limit. Thus, if we consider a gauge parameter ^ in the contin- 
uum, the lattice quantity a^gQA^{x) is generated from a Gaussian distribution with width 
o = ^/2N^^Jp, instead of a width ^/^. 

Note that a = -y/^ if j8 = 2Nc and that for j8 < 2Nc the lattice width a is larger 
than the continuum width ^/^, making the convergence problem discussed above more 
severe. Thus, in the SU(2) case, one has o = only for j8 = 4, corresponding to a 
lattice spacing a ^ 0.001 fm. On the contrary, in the SU(3) case, one has a = \fh, for 
j8 = 6, corresponding to a = 0. 102 fm. Also, for a fixed t'Hooft coupling ^Nc, we have 



TABLE 3. Values of j3 for which the numerical 
gauge- fixing algorithm showed convergence. Re- 
sults are reported for three different gauge groups 
and four different lattice volumes. In all cases the 
gauge parameter ^ was 1 (Feynman gauge). *[In 
these two cases only a few configurations have 
been considered and more tests are needed.] 
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j8 oc AT^ and o <^ ^/TJn^.. This suggests that simulations for the linear covariant gauge 
are probably easier in the SU(A^t ) case for large Nc- 

In order to test this hypothesis we simulated the SU(2), SU(3) and SU(4) cases for 
a gauge parameter ^ = I and lattice volumes V = 8"*, 16'*, 24"^, 32"^ for the values of j8 
reported in Table 2. They correspond, respectively, to a t'Hooft coupling g^Nc = 8/3 
(j8i), to a plaquette average value of about 0.65 (^2) and of about 0.6 (f^s) and to a 
string tension (in lattice units) of about a^o — 0.044 (^84), giving a ^ 0.09 fm. In Table 
3 we present, for each pair (S\J (Nc), V), the values of the lattice coupling j8 for which 
the gauge-fixing algorithm showed a numerical convergence.^ One clearly sees that the 
situation improves when the number of colors Nc is larger. 

CONCLUSIONS 

We have recently introduced a minimizing functional for the linear covariant gauge 
which is a simple generalization of the Landau-gauge functional. The new approach 

solves most problems encountered in earlier implementations and ensures a good quality 
for the gauge fixing. Here we have shown that, by using different discretizations for the 
gluon field (such as the angle projection) and by considering a gauge group SU(A^c) with 
Nc sufficiently large — i.e. A'c = 4 or maybe even Nc = 3 — one should be able to do 
simulations for ^ = I (Feynman gauge) and for large lattice volumes (in physical units). 
Let us note that a numerical study of the infrared behavior of propagators and vertices 
at^^O could provide important inputs for analytic studies based on Dyson-Schwinger 
equations [28, 29]. Moreover, it has been proven [30, 31, 32] that the background-field 
Feynman gauge is equivalent (to all orders in perturbation theory) to the pinch technique 
[32, 33]. Thus, numerical studies using the Feynman gauge, which corresponds to the 
value t, — I, could allow a nonperturbative evaluation of the gauge-invariant off-shell 
Green functions of the pinch technique [34]. 



For these tests we used the standard overrelaxation algorithm [22, 23, 24]. 
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